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ABSTRACT 

We have detected asymmetry in the symbiotic star CH Cyg through the measurement of preci- 
sion closure-phase with the IONIC beam combiner, at the IOTA interferometer. The position 
of the asymmetry changes with time and is correlated with the phase of the 2.1-yr period 
found in the radial velocity measurements for this star. We can model the time-dependent 
asymmetry either as the orbit of a low-mass companion around the M giant or as an asym- 
metric, 20% change in brightness across the M giant. We do not detect a change in the size 
of the star during a 3 year monitoring period neither with respect to time nor with respect 
to wavelength. We find a spherical dust-shell with an emission size of 2.2±0.1 FWHM 
around the M giant star. The star to dust flux ratio is estimated to be 11.63±0.3. While the 
most likely explanation for the 20% change in brightness is non-radial pulsation we argue 
that a low-mass companion in close orbit could be the physical cause of the pulsation. The 
combined effect of pulsation and low-mass companion could explain the behaviour revealed 
by the radial-velocity curves and the time-dependent asymmetry detected in the closure- 
phase data. If CH Cyg is a typical long secondary period variable then these variations could 
be explained by the effect of an orbiting low-mass companion on the primary star. 

Key words: binaries: symbiotic - stars: imaging - stars: individual: CH Cygni - techniques: 
high angular resolution - techniques: interferometric. 



1 INTRODUCTION 

Symbiotic stars are objects presenting combination spectra of a hot 
ionised nebula and the cool continuum absorption molecular fea- 
tures of a late-type star. Nowadays, symbiotic stars are understood 
as mass-transfer binaries of short period, from a few to 10 years. 
The separation can vary from a few AU to slightly more than 10 
AU. The symbiotic pair is usually composed of a cool giant star 
with an accreting compact object, either a white dwarf or a neutron 
star. 

CH Cyg is one of the most studied of symbiotic variables. 



* Affiliated to Scottish universities physics alliance (SUPA). 
© 2008 RAS 



The star presents a composite spectrum of a M6-7 giant star during 
quiescent phase and a hot component blue continuum from 6000 
to 9000 K temperature and low exci tation line spec trum during the 
active phase foeutsch et al.lll974) . IWebster & AllenI ^9751) classi- 
fied the star as an S-type symbiotic with no hot dust, but long term 
multi-wavelength photometry study of the star jTaranova & ludinl 
119881) has shown that hot dust appeared in the system after the 1984 
outburst. Th e dust was modelled as a sphe rical shell of inner radius 
of 15 AU bv lBogdanov & Taranoval ( 1200 ll) through spectral energy 
distr ibution (SEP) fitti ng. 

iDvck et al] j 19981) measured an angular diameter of 10.4 mas 
at 2.2 /^m for CH Cyg with infrared interferometry. lYoung et al.l 
I2OOO) observed CH Cyg with the Cambridge optical aperture syn- 
thesis telescope (COAST) in 1999. The obtained visibility and 
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Table 1. Log of observations. The IOTA 3T telescope configuration refers to the location of the A, B, C telescopes along the NE, SE and NE arms. 



Date Mean Phase Telescope A AA R''^ Calibrator 

(UT) JD (p™) (^'^™) names 



2004Apr23 


2453119 


0.40 


IOTA 


A35- 


-B15- 


-ClO 


1.51 


0.090 




a Lyr, a Aql 


2004Apr24 


2453120 


0.40 


IOTA 


A35- 


-B15- 


-CIO 


1.64 


0.100 




a Lyr, a Aql 


2004Apr25 


2453121 


0.40 


IOTA 


A35- 


-B15- 


-CIO 


1.64 


0.100 




a Lyr, a Aql 


2004Apr26 


2453122 


0.40 


IOTA 


A35- 


-B15- 


-CIO 


1.78 


0.090 




a Lyr, a Aql 


2004Apr29 


2453125 


0.41 


IOTA 


A35- 


-B15- 


-CIO 


1.78 


0.090 




a Lyr, p Ser 


2004Apr30 


2453126 


0.41 


IOTA 


A35- 


-B15- 


-CIO 


1.78 


0.090 




a Lyr, u Hya 


2004May01 


2453127 


0.41 


IOTA 


A35- 


-B15 


-CIO 


1.78 


0.090 




a Lyr, a Aql 


2004Sep04 


2453105 


0.38 


Keck A, Golay mask 


1.64 


0.025 




a Lyr 


2005Jun06 


2453528 


0.94 


IOTA 


A35- 


-B15- 


-CIO 


1.66 


0.300 




a Lyr 


2005Jun08 


2453530 


0.95 


IOTA 


A25- 


-B15- 


-CIO 


1.66 


0.300 




a Lyr 


2006Apr24 


2453850 


0.37 


IOTA 


A35- 


-B15- 


-CIO 


1.66 


0.300 


39 


a Lyr, /3Her 


2006Apr30 


2453856 


0.38 


IOTA 


A35- 


-B15- 


-CIO 


1.66 


0.300 


39 


a Lyr, /3Her 


2006May01 


2453857 


0.38 


IOTA 


A35- 


-B15- 


-CIO 


1.66 


0.300 


39 


a Lyr, /3Her 


2006May02 


2453858 


0.38 


IOTA 


A35- 


-B15- 


-CIO 


1.66 


0.300 


39 


a Lyr, /3Her 



°' Only applicable to the IOTA spectrograph. 



closure-phase data were best modelled by an elliptical, limb- 
darkened star. Their interpretation of these findings was that the 
ellipticity of the star was either due to an extension of the M giant 
atmosphere or to the partial eclipse of an orbiting red giant com- 
panion as proposed bv lSkopal et al] ( ll996h . 

CH Cyg shows photometric and radial-velocity variations. 
In photometry two main periods ar e found: a small ampli- 
tude (0.1 magnitude) '^lOO day period jMikolaiewski et aT 
likely caused by stellar pulsation and a ~770 day period 
(Mikolajewski et al. 1992) or a ~1 magnitude, ~750 day period 
dSkopal e t al. 2007). T hese p hotometric variations are not always 
detected (iMunari et al] 199q) and ar e not related to the 100-day 
pulsation period. iHinkle et al] ( Il993h pointed out that the photo- 
metric variations of CH Cyg are far longer than the fundamen- 
tal pulsation mode for th is star, which is a first-overtone pulsator 
jMikolaiew ski et al. 1992). The radial-velocity variations from the 
literature tHinkle et aL 1993i) show two periods: a 15.6-yr long pe- 
riod and a 2.1-yr (750 days) short period. 

iHinkle et al.l j 19931 1 20081) suggested a correlation between 
the 750/770 days photometric period and the 750 days radial- 
velocity period. They also remarked that the photometric varia- 
tions of CH Cyg are simil ar to those found in long secondary pe- 
riod (LSP) variable stars ( IHinkle et al.ll2006h . This type of vari- 
ability is found in some semi-regular variables and in about the 
25% of the L arge Magellani c Cloud (LMC) se mi-regular vari- 
ables. [ginSeet al. (2002) and lWood et"ai] ( l2004t) found that sev- 
eral semi-regular variables also show spectroscopic behaviour con- 
sistent with LSP variability. The cause of the LSP is currently 
unkno wn but possible explanations are highlighted in I Wood et al.l 
( |2004|) . They conclude that the most likely explanation for LSP 
variations is a low order non radial pulsation on the outer radiative 
layers of the giant star. 

I Wood et"aLl ( l 19991) also discovered that LSP variables follow 
a period-lumino s ity (P -L) relation which he called "sequence D". 
ISoszvnski et al.l ( l2004b noted that Wood sequence-D variables 
overlap with sequence-E contact binaries, implying that sequence- 
D is indeed a class of binaries. ,S^zyhski ( 2007) also found that 
5% of LSPs in the LMC present ellipsoidal-like or eclipsing-like 



modulation that are usually shifted in phase with respect to LSP 
light curves. 



IHinkle et al. I h993l) proposed a model where CH Cyg was a 
triple system with the symbiotic pair in a 2.1-yr orbit. The rea- 
sons for having the symbiotic pair on the 2.1-yr orbit were that 
no known S-type symbiotic star had orbital period larger than 5 
years, the 2. 1-yr period was too long for a M giant fundamental- 
mode pulsation and there was weak evidence for a high inclination 
15.6-yr orbit. The third star was e ither regarded as a G-K dwarf 
( lHinkleetalJll993h or a M giant dSkopal et al]|l996l) . The incli- 
nation of the 15.6-yr orbit was unknown at the time but was re- 
cently infe rred from the several eclipses reporte d in the literature 
dMiko l aiewski et al.lll 987: Eyres et al. 2002; S okoloski & KenvonI 
1 20031) . ISchmidt et alT (2006) suggested that the 2.1-yr period was 
caused by a pulsation in the M giant and not by a close binary. 

Ther e is controversy on the shape of the possible orbit of the 
close pair. lHinkle et al. 1(11993') argued that the asymmetric line pro- 
files could be caused by a M giant star irradiated by a white dwarf. 
An asymmetric line profile could lead to a false elliptic solu tion for 
an orb it obtained from radial velocities. According to Hinkle et al.l 
( Il993t) the orbit of CH Cyg should be circular due to tidal interac- 
tion with the M giant. 



IHinkle et al.l j2008l) re-examined the conclusions of the 
IHinkle et al. I h993h paper. They concluded that the 2.1-yr veloc- 
ity variation is consistent with LSP variation and that the white 
dwarf responsible for the activity in the system is on the 15.6-yr 
orbit. The 2. 1-yr period would be caused either by non-radial pul- 
sation of the star or by a low-mass companion in close orbit to the 
M giant. 

This paper presents the results of infrared interferometric ob- 
servations per formed in 2004-20 06 at the infrared optical telescope 
array (IOTA) dTraub et al.ll2004 and at the Keck-1 telescope fitted 
with an aperture mask. The main aim of this paper is to provide 
unique observational data that could help to understand the nature 
of the mysterious 2. 1-yr oscillation in radial velocity for this star. 
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Table 2. Calibrator information. 



Calibrator 


Spectral 


Adopted UD 


Reference(s) 


name 


type 


(mas^ 




a Lyr 


AOV 


3.22±0.01 


Absil et al. (2006) 


/3Her 


G7IIIa 


3.40±0.03 


This work 


p Ser 


K5III 


3.28±0.04 


Borde et al. (2002) 


a Aql 


A7V 


3.46±0.04 


van Belle et al. (2001) 



milliarcseconds 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were performed at the IOTA interferometer and at 
the Keck-1 telescope. IOTA was a long-baseline optical inter- 
ferometer located at the Smithsonian Institution's Whipple Ob- 
servatory on Mount Hopkins, AZ. IOTA operated from 1995- 
2006, an d was used as a t e stbed for new cuttin g-edge tech- 
nologies teerger et all 1200 ll : iMonnier et al] |2003|) . IOTA pro- 
duced a l arge number of astronomy results over the past few 
years jMennesson et al. 2002; Ohnaka et al. 2003; Monnier et al.l 



2004; Perrinetal. 2004; Kraus et al. 2005; Millan-Gabet et al. 
2006.; Monnier et al... 2006; Zhao et al.. .2007. ; .Kraus et aL ,2007. ; 



Ragland et al.ll2008l ; lLacour et alj|2008h . 



Observations performed at the Keck-1 telescope used the near 
infrared camera (NIRC) and an aperture mask that converted the 
telescope pupil into a sparse interferometric array of 9.8 m maxi- 
mum baseline. For detailed discussio n of the Keck apertu re-mask 
experiment and scientific rationale see iTuthill etal.l ( l2000l) . 

In Table [T] we present a journal of our observations, listing 
date, filter and calibrator star. In Table[2]we detail the physical prop- 
erties of the calibrat ors. Our H-band data used the IONIC combiner 
teerger et al]|2003h with narrow H-band filters at the IOTA inter- 
ferometer and at the Keck telescope for the observations of 2004, 
while data was acquired using a standard H-band filter at IOTA for 
2005. Data from 2006 used a low-dispersion spectrograph which 
provided seven spectral channels across the H-band wi t h an R =39. 
For a descript i on of the spectrograph see lRaglandetal] ( l2003h and 
IPedretti et al.l (|2008|). First re sults with the spectrograph were pub- 
lished bv lLacour et all fcoOSh . 

For th e aperture-maski ng e xperiment we refer the readers to 
the work of lMonnieJ ( 1 1 9991) and lTuthill et~^ feOOd) for the data- 
analysis, procedures adopted to extract visibilities and closure- 
phases in OIFITS format (Pauls et al. 2005). 

The data reduction pipeline f or the IONIC combiner was de- 
scribed in detail in .Monnier et'ai] j2004l) . Briefly, reduction of the 
squared visibilities (V^) follow ed the same method explained by 
ICoude Du Foresto et al.l l ll997h . Interferograms were corrected for 
intensity fluctuations and bias terms from readout noise and photon 
noise. The power spectrum of each interferogram was calculated in 
order to measure V'. A transfer matrix was used to take in account 
the variable flux ratio f or each baseline. The absolute calibration 
accuracy was studied bv lMonnieretal.l ( l2004l) by observing single 
stars of known size. 

Closure phases for the IONIC combiner were obtained using 
two in dependent methods; one was developed by [Baldwin et al.l 
l l 19961) for the Cambrid ge optical apert ure synthesis array 
(COAST), and the other by HaleetalJ ( |2003|) for the infrared spa- 
tial interferometer (ISI). In order to measure meaningful closure 
phase, fringes must at least be present in three baselines and the 
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Figure 1. The uv coverage for CH Cyg for 2004 (circles) 2005 (squares) 
and 2006 (triangles). The shades of grey represent the wavelength of the 
data points. 



fringe packets must overlap, to be detected in the same coherence 
time. The largest error in closure phase offset for a point source 
was caused by chromaticity in the combiner which limits the abso- 
lute precision when source and calib rator are not the same spectral 
type. Engineering tests performed bvl Monnier et al.l ( l2004h showed 
that the closure phase varied systematically by 1.4 ± 0.3° between 
a cool star of spectral type M3 and a hot B8 star. 

The IOTA data pipeline produced visibility and closure phases 
in OIFITS format, which can be easily imported in imaging or mod- 
elling programmes using libraries provided by John Young, for C 
and PythorQand John Monnier for IDL0. A standard 2% systematic 
error was added in quadrature to the visibility and closure-phase 
data as in lMonnier et al.l j2004h . Calibrated data in OIFITS format 
will be made available on request for interested investigators. 



3 MODELLING 

IWoode7ZI ( l2004 l) conducted a thorough review on the causes 
of the LSP variations in asymptotic giant branch (AGB) stars. They 
ruled out several possible models, among which were radial pul- 
sation, companion in close orbit, spots on the star and modulation 
from an ellipsoidal-shaped AGB star. They concluded that non- 
radial pulsation was the most likely explanation for LSP. In their 
recent paper Hinkle et al. (2008) after a thorough review of the lit- 
erature on CH Cyg applied a similar approach to rule out possible 
models explaining the 2. 1-yr change in radial velocity in CH Cyg. 
We used our interferometry data to verify some of the hypotheses 
discussed in these papers. Due to limited uv -plane coverage of the 
data (see Figure [TJ, in particular for the 2005 epoch, we could not 
resort to model-independent imaging of the CH Cyg system. For 
this reason we used parametric modelling to derive the size of the 
star, the FWHM size of the dust and the position and distance of the 



^ http://www.mrao.cam.ac.uk/research/OAS/oi_data/oifits.html 
^ http://www.astro.lsa.umich.edu/ monnier/oi.data/index.html 
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Figure 2. Visibility plot comprising all data from 2004 to 2006. A simple 
uniform disc plus a Gaussian disc model representing the star and the dust 
are superposed to the data. The contributions from the dust (dashed line) 
and the star (dash-dotted line) are also over-plotted for clarity. 

asymmetries detected in the closure-phase data. For model fitting 
we used publicly available least-squ ares minimisation rout ine jf[ 

Our modelling was similar to iRagland et all (l2008h except 
that we did not need to model multi-wavelength sizes for the star, 
since CH Cyg does not change size appreciably with respect to 
wavelength. We decided to test the following hypothesis to inves- 
tigate the cause of the LSP variations and to interpret our data: (1) 
radial pulsation of the star, (2) presence of dust inside the 15.6-yr 
orbit (3) spots on the star, (4) M giant companion in a 15.6~yr or- 
bit, (5) dwarf companion in a 2. 1-yr orbit, (6) non-radial pulsation. 

3.1 Radial pulsation and dust 

In order to test (1) and (2) a simple model composed of a uniform 
disc (UD) for the star and a Gaussian disc (GD) for the dust was 
first attempted in order to obtain a size for the star and for the dust. 
All data from all epochs were used for this model since, by visual 
inspection, our visibility points superposed quite well, indicating 
that the size of the star did not change appreciably outside the er- 
ror bars of the data, neither with time nor wavelength nor position 
angle. 

Figure [2] shows the result of the fit. The data were smoothed 
using an azimuthal average due to the otherwise very large num- 
ber of data points present on the graph. For each bin we used the 
mean of the original data points weighted by their errors. The er- 
ror on each new data point was the standard deviation for the bin. 
The fit was performed on the original and non-smoothed data. 
Table |3] shows the parameters obtained from the fit. The value of 
8.74±0.02 milliarcseconds (from now on mas), for the diameter of 
the M giant is the most accurate so far thanks to the large amount of 
data used. This value is close to the valu e of 7.8±0.6 obtained with 
infrared interferometry in June 2001 bv lHofmann et alj ( |2003h us- 
ing a simple UD fit. The errors were derived using bootstrap statis- 
tics on the data set. The full-width half maximum size of the Gaus- 
sian dust emission was 19.13±1.00 mas or 2.2±0.1 stellar diam- 
eters FWHM, showing that hot dust exists close to the M giant. 

A non-linear least squares curve fitting (MPFIT), developed by Craig 
Marquardt. http://cow.physics.wisc.edu/craigm/idiyi 



Table 3. Size of star and size of the dust-shell emission. We used a uniform 
disc to model the star and a Gaussian disc to model the dust. A model with a 
spherical dust-emission and a model with an elUptical dust-emission were 
attempted. The two models produced a very similar reduced 



Size 


FWHM 


Flux ratio 


P.A. 


Axis ratio 


Reduced 


M giant 


dust_^ 


(M giant 




dust 




(mas) 


(mas) 


/dust) 


(°) 


(M/m) 




8.74±0.02 


19.1±1.0 


11.6±0.3 


0.0 


1.0 


1.3 


8.74±0.02 


19.2±0.9 


11.6±0.1 


103±5 


1.28±0.01 


1.2 



Size of the dust emission. 

A marginally improved reduced was obtained by fitting an el- 
liptical dust distribution around the star. However the difference in 
reduced was too small in order to justify an asymmetric model 
for the dust emission in the near infrared. The parameters from the 
ellip tical dust-emission mod el are also listed in Table|3] 

[Thompson et al.l ( l2002h monitored an oxygen-rich and a 
carbon-rich Mira star measuring the change of angular size with 
respect to the pulsation cycle at the Palomar testbed interferometer 
(PTI). We did not detect any such change in CH Cyg. Unfortunately 
our coverage of the 2.1-yr period was quite limited (basically two 
points at phase 0.4 and one point at phase 0.9 of the "orbital" pe- 
riod). This coverage is insufficient to comple tely rule out rad i al pul- 
sation for this star. However, we notice that lHofmann et al] ( l2003h 
obtained a diameter of 7.8±0.6 mas in June 2001, using a simple 
UD model with three visibility points. Considering the crude UD 
model used that does not take in consideration the dust shell, this 
diameter is not very different from our measurement and would 
indicate that the star did not c hange diameter with time. Also, ra- 
dial pulsation was ruled-out bv lHinkleetalJ ( l2008h as the cause for 
secondary peri od in CH Cyg since the period-luminosity relation 
for AGB stars faughes & Wooj|l990l) would produce a period of 
about 250 days for a K=-7.5 star not 770 days. 

3.2 Spots on the star 

In order to model the closure-phase signal expected from a spot- 
ted star, we used an additional uniform disc that could be placed 
at different position angles and separation from the centre of the 
UD -I- GD model representing the M giant and the dust emission. 
Flux ratios between the M giant and the companion/spot and be- 
tween the M giant and the dust were allowed to change. The size 
of the additional UD was also free to change. The three epochs 
were analysed separately in order to detect asymmetries that would 
change with the observing epoch. We performed a parameter-space 
search in an attempt to identify the position of the asymmetry. 

We could not find any solution with an unresolved or mod- 
erately resolved spot on the surface of the star. All the solutions 
converged to structure outside the disc of the M giant unless we re- 
stricted the flux to 10% or more of the flux of the M giant as done 
in Section l331 

3.3 M giant companion on the 15.6-yr orbit 

We tested the hypothesis that the companion is a red gia nt in the 
15.6- yr orbit as discussed in the model proposed by Skopal et al.l 
( Il996h . That model was devised to explain the eclipses observed 

© 2008 RAS, MNRAS 000,[T]-?? 
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Table 4. Orbital positions for the low-mass object. 



Mean 
JD 



Sep 
(mas) 



P.A 

(°) 



Flux ratio 
M giant/ 
dwaii' 



Flux ratio 
M giant/ 
dust 



Reduced 



2453122.9 


7 +3 

' -2 


188 


78±1 


9.6±0.4 


0.9 


2453529.0 


8+3 
= -2 


356 


+ 4 
-4 


88±5 


18.7±1.0 


0.2 


2453529.0 


32tt 


331 


f3 
-3 


104±4 


18.4±0.8 


0.3 


2453855.2 




211+37 
^'^-42 


74±1 


14.2±0.4 


0.9 



in the 15.6-yr orbit and kept the symbiotic pair in the 2. 1-yr or- 
bit, given that there are no known symbiotic stars found on an 
orbit of period as long as 15.6 yrs. We restricted the flux ratio 
of the M giant/companion to values around 8.6, as expected in 
iTaranova & Shenavrin (200J). The field of view (FOV) of IOTA 
was limited by the bandwidth of the photometric filter used: 
FOV = A^/AAB, where A is the wavelength AA the bandwidth 
and B the baseline. For the largest bandwidth used (0.3 /im at 1.65 
^m) and a baseline of 38 metres we obtained a minimum FOV of 
50 mas. We performed a 50 mas wide search in all our data sets. 

We did not find any trace of a companion in our best data sets 
of 2004 and 2006 when using the Taranova flux ratio. A second 
red giant should have been evident in the data. In particular the 
Keck telescope aperture-masking experiment should have easily 
detected a sec ond giant star down to a flux ratio M gi ant/companion 
of about 100 jlreland et alj2008l : lKraus et alj2008l) . 

3.4 Dwarf companion in a 2.1-yr orbit 

We tested the hypothesis of a faint companion orbiting the M giant 
as in Section [T2I and Section [33] We restricted the flux contribu- 
tion of the companion to less than 2% of the total flux in order 
to simulate a large Am between the M giant and the companion. 
As a consequence we found asymmetries outside the M star in all 
three data sets. Figure [3] shows the likelihood maps obtained from 
the reduced surfaces. The dotted-line ellipses are the errors on 
the positions of the companion which are quite large in the East- 
West direction due to the limited uv coverage of the IOTA in that 
direction. 

Table |4] lists the separations, position angles, flux ratios and 
reduced of the asymmetries for the three epochs. The UD size 
of the companion converged to a point source for all epochs. The 
error bars on the parameters were derived from the error ellipses. 
The 2005 data set converged to two separate solutions: one at 8 mas 
separation and another 32 mas separation. The second solution was 
likely due to a degeneracy caused by the limited amount of data 
available for the 2005 epoch. We could fit the 2. 1-yr elliptical orbit 
to the 32 mas position, with a of 0.3. The semi-major axis of this 
orbit was 25.6±0.8 mas which produced a far smaller luminosity 
and a much shorter distance than expected for this star. For this rea- 
son w e excluded this solution. We must point out that lBalega et al.l 
( l2007h detected a faint companion with speckle interferometry in 
2004 at43±l mas separation and 24. 1° ±2.1 "position angle. How- 
ever we do not believe that our 32 mas position is related to this 
detection. In fact the Keck telescope aperture-masking experiment 
should have easily detected a companion down to a flux ratio M gi- 
ant/companion of about 100 in our 2004 data jlreland et al]|2008l : 
iKraus et al.ll2008h but such detection did not happen. 

In order to investigate the hypotheses that the detected asym- 
metries were the signature of a faint companion in a 2. 1-yr orbit 



we attempted orbit fits to the astrometric positions derived from the 
IOTA closur e-phase data using infrared radial-velocity orbital so- 
lutions from [llinkle et al.l ll2 0081. We attempted orbi tal fits using 
a circular orbit I Fekell 20081) which was discarded in iHinkle et alj 
I2OO8) due to the large resid uals in the orbit fit. W e also used the 
elliptical orbit solution from H inkle et al.l ( l2008h . We obtained a 
reduced of 0.1 for the circular orbit and a reduced of 0.3 
for the elliptical orbit. Such small reduced values are possible 
given the small number of degrees of freedom (3 from the 6 data 
points and 3 free parameters) and the quite large error bars on the 
astrometric positions. The obtained orbits are shown in Figure |4] 
The combined orbital parameters from radial velocity and interfer- 
ometry and some derived parameters are shown in Table [5] The 
first part of the table lists the orbital parameters from lHinkle et al.l 
([2008) and Fekel (2008). The errors on the parameters obtained 
from interferometry were derived using Monte Carlo simulations. 



3.5 Non-radial pulsation 

Although we could not obtain acceptable reduced x^ from the 
"spots on the star" model we managed to obtain reasonably good 
fits in simulations of large asymmetries on the star. We used low- 
order spherical harmonics to simulate large flux variations (up to 
20%) across the star. Such a dramatic brightness change could sim- 
ulate the closure phase signal of CH Cyg but seemed an unlikely ex- 
planation even in term of non radial pulsation: we are not aware of 
any physical mechanism that could produce such a dramatic change 
of brightness across a star. Figure [5] shows the models and the cor- 
responding fits of visibility and closure phase. The reduced was 
reasonably close to the dwarf-companion model. The asymmetry 
appears to rotate with the 2. 1-yr period. 



4 DISCUSSION 

iHinkle et al.l ( |2008|) restricted the possible explanation for the 2. 1- 
yr secondary period of CH Cyg to a low-order g-mode non-radial 
pulsation of the M giant or a to low-mass companion (0.2 Mq ) in 
close orbit to the M giant. In this model the companion responsible 
for the activity is on the I5.6-y r orbit. 

According to lHinkle et al. I1I2OO8) a 0.2 Mq companion would 
have a temperature of about 3200K and would be spectroscopically 
indistinguishable from the M giant. Since CH Cyg is single-lined 
binary/triple star the masses of the components cannot be derived 
directly. We can however test the derived parameters against the 
published literature, assuming the mass of one of the components. 
Table [6] shows the change of the derived parameters for different 
values of the mass of the companion. Mi is the mass of the M giant 
in solar masses, Ri the radius in solar radii, Li the luminosity in 
solar luminosities. The semi-major axis "a" of the orbit in physical 
units of AU was obtained using Kepler's law. D is the distance in 
parsecs. The table is divided in two parts, one concerned with the 
circular orbital solution and o ne with the el l iptica l solution. 

The mass of 2 M(7) fromlHinkle et alj ( |2008|) . a luminosity of 
6900 L0 from Biller et al. ('200^ and a radius of 280±65 R© ob- 
tained by Schild et al. ( 1999) through infrared spectroscopy were 
used to restrict the solutions listed in Table |6] A value of 
0.32 M0 for the low-mass companion yielded a mass of 2 Mq a 
radius of 250 R© and a luminosity of 6517 Lq for the M giant, very 
close to the values from the literature. 

Also, the distance obtained from the size of the circular or- 
bit derived using Kepler's law and the apparent size of the orbit 
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Figure 3. Binary-star search in space. A simple binaiy-star model was fitted to the visibility and closure-phase data from different epochs. The M giant 
was kept at the centre of the field while the companion was placed in all possible positions of a 40 X 40 mas square grid. A reduced value was obtained for 
each position and the values recorded on a two-dimensional array. Left column shows a likelihood surfaces derived from the arrays for our 2004, 2005a and 
2005b (non-unique solutions) and 2006 data. The positions of the asymmetries are encoded in the likelihood map. The white crosses represent the positions 
of the asymmetry and the white ellipses encode the uncertainty of the position. The centre column shows data-versus-model plots for closure phase. The right 
column shows visibihty-versus-spatial-frequency plots (filled circles) with superposed points derived from the model (open diamonds). Note that the high 
density of data points makes it very hard to distinguish between filled circles and open diamonds. We observed that the closure phase flipped sign between 
2004 and 2005 and between 2005 and 2006 meaning that the the detected asymmetry was in the opposite direction in 2005. 
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Figure 4. The astrometric orbit of CH Cyg. The plots show an elhptical orbit fit (left) and a circulai' orbit fit (right). Superposed to the orbits are a uniform 
disc representing the star and a Gaussian disc representing the dust emission in the system. The flux contribution from the dust was exaggerated to render the 
dust extent visible in the picture. The diamonds are the expected positions of the companion relative to the M giant, according to the ephemeris. The observed 
positions of the secondary component are marked with error ellipses (dotted line) centred around a star symbol. The triangles are the expected positions of the 
companion during the observations at the COAST inteii'erometer. The dashed-line ellipses represent an elliptical model of CH Cyg from Young et al. (2000). 
The major axis of the ellipses is also shown to better appreciate the orientation of the ellipses. 



Table 5. Orbital parameters of the possible low-mass object. 



Parameters 


Circular 


Elliptical 




solution 


solution 




Radial velocity 




P (days) 


749.8±2.3 


750.1±1.3 


TO (HJD) 


2446823.2±7.7 


2447293.5±12.9 


Lun 


0.0 


229.5±7.7 


e 


0.0 


0.330±0.041 


K(Km s-i) 


2.87±0.13 


2.87±0.13 


7 (Km s-i) 


-59.93±0.10 


-59.91±0.09 


a sin I (Km) 


2.96xl0'±0.29xl0^ 


2.79xl0'±1.23xl0' 


f(m) (Mq ) 


0.00018±0.0002 


0.00015±0.0002 




Interferometry 




in 


138±10 


146±6 


nn 


347±7 


337±8 


a (mas) 


7.1±0.3 


6.3±0.3 



in milliarcseconds, was 296 pc, comparable within errors to the 
244^35 pc distance ob tained from the revi sited data reduction of 



the Hipparcos paralla x i van Leeuweiill2007h . 

ISoszvnski et al. I j2004h found that Wood's sequence-D vari- 
ables are a co n tinuat ion of sequence-E ellipsoidal variables. 
ISoszvnski et al.l ( l2007h require that the ratio between the ra- 
dius of the star and the semi-major axis of the orbit should be 
R/a ~ 0.4 for the binary explanation of the LSP (Equation 5 in 



ISoszvnski et al.l l2007*)). For our circular orbit solution the ratio de- 
rived by the angular diameter of the star and the semi-major axis of 
the hypothetical orbit is 0.6, ~50% larger than the required value. 

Soszvhski ( 2007) also proposed a model where the LSP vari- 
ation are caused by mass loss from the giant to the low-mass com- 
panion. Since we detected hot-dust emission inside the possible 
2.1-yr orbit we cannot exclude that LSP photometric variation are 
caused by dust trailing the low-mass companio n. There have been 
claims of ec lipses in the 2. 1-yr orbit of CH Cyg ( ISkoDaletalJl996l : 
lliiimal 19981) . The inclination of the circular orbit obtained from the 
interferometric data would prevent eclipses but if the dust is clumpy 
and is trailing the companion it could be responsible for occasional 
photometric variations and could simulate eclipses. 

Another clue in favour of the low-mass companion expla- 
nation for the LSP variation of CH Cyg comes from the inde- 
pendent interferometric observations of Young et al. (2000) at the 
wavelength of 905 nm. The visibility and closure-phase data from 
the COAST interferometer were modelled by an elliptical, limb- 
darkened star (parameters from that model are reproduced in Ta- 
ble |7). Figure |4] shows the astrometric orbit of CH Cyg for the el- 
liptical and circular orbit solution superposed to the models from 
lYoung etal.l l l2000 ). Interestingly the minor axis of the two ellipses 
is very close to the radius of the M giant obtained from our model, 
while the major axis intersects the predicted orbital position of the 
companion on the circular orbit. The major axis of the ellipses did 
not intersect the orbital positions of the companion on the elliptical 
orbit. Also, our elliptical orbit fit to our astrometric positions had a 
3-times worse than the circular orbit fit. 

According to lHinkle et alj ( |2008|) the most powerful argument 
against the close-binary explanation of LSP variation in CH Cyg is 
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Figure 5. The 2004 to 2006 visibility and closure-phase data was also fit by an asymmetric brightness distribution on the surface of the star. The left column 
shows models of the star for different epochs. The asymmetric flux distribution was modelled using spherical harmonics. The centre column shows data-versus- 
model plots for the closure phase. The light column shows visibility-versus-spatial-frequency plots (filled circles) wi th superposed points derived from th e 
model (open diamonds). The asymmetric change in biightness across the star could be caused by a non- radial pulsation JWood et alj2004l : lHinkle et alj2008l) . 
However we are not aware of any pulsation mechanism that would produce a 20% change in brightness across the star. The closure phase changed sign and 
the asymmetry flipped of 180° in 2005. 



the shape of the radial-velocity curve: the most likely orbit for the 
low-mass companion would be elliptical due to the large radial- 
velocity residuals obtained from fitting a circular orbit. On the 
other hand Hinkle et al. (2008) derived a mass of 2.0 for the 
M star, based on evolutionary arguments and argued that the Roche 
lobe for a 2 Mq - 0.2 Mq binary would constantly change from 
(1 + e)a at apoastron to (1 — e)a at periastron for an elliptical 
orbit. A 280 R© giant would fill the Roche lobe at each perias- 
tron passage generating large mass loss. We argue that distortion 



by proximity effect teatoij|2008l) could also change the shape of 
the radial velocity curve and therefore the circular orbit solution 
cannot be eliminated on the basis of this argument. 

iHinkle et al. I( l2008h argue that non-radial pulsation could also 
reproduce the observed radial velocities and photometric variation 
of CH Cyg. The problem with the non radial pulsation argument is 
that low-order g modes are evanescent in convective regions and 
there is no known physical mechanism that could explain the non- 
radial pulsation for M giant stars where radiative transfer is mostly 
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Table 6. Derived red giant parameters as a function of the faint compan- 
ion's mass. 



Ma (Mq ) 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 






Circular orbit 








Ml (Mq ) 


0.3 


0.9 


1.8 


2.8 


3.9 


5.2 




156 


221 


270 


312 


349 


382 


Li (Lo ) 


2037 


4073 


6110 


8146 


10183 


12219 


a(AU) 


1.2 


1.7 


2.0 


2.3 


2.6 


2.9 


D(pc) 


166 


235 


288 


332 


371 


407 






Elliptical orbit 








Ml (Mq ) 


0.2 


0.7 


1.4 


2.2 


3.2 


4.2 


Rl (R© ) 


165 


233 


285 


330 


369 


404 


Li (Lq ) 


2267 


4534 


6801 


9068 


11335 


13601 


a(AU) 


1.5 


1.6 


1.9 


2.2 


2.5 


2.7 


D(pc) 


143 


248 


304 


351 


392 


430 



Table 7. Elliptical star model parameters from Young et al. (2000). 



Epoch 


Major axis 


Axial ratio 


P.A 




(mas) 




(°) 


99/08 


11.5 ±0.2 


0.84 ±0.05 


126±9 


99/09 


11.2 ±0.2 


0.79 ±0.03 


136±5 



convective. As we show in Figure [5] an asymmetric brightness dis- 
tribution on the surface of the star could also reproduce our ob- 
served closure-phase signature. The flux variation across the star 
must be very large (20%) in order to explain the closure-phase re- 
sults and we are not aware of any physical mechanism that could 
produce such a dramatic change of brightness across a star. 

Close encounter with another object can produce n on-radial 
oscil lations on a fluid star through tides according to lEriguchil 
( Il990 ). Circularisation of early-type main sequence binaries are 
also Imown to cause non-radial g-mode resonant oscillations. In 
late type stars turbulence is very effi cient in damping these oscil- 
lation and circularising binary orbits. Ilvanov & PapaloizoJ j2004l) 
studied tides in fully convective stars. They came to the conclusion 
that resonant tides may be possible in fully convective stars. 

CH Cyg is a very complex object. If it is a triple system the 
interactions of the companions with the M star could be very com- 
plex. The signature in the radial velocity could be caused by a com- 
bination of movement of the low-mass companion and non-radial 
pulsation of the star. This could explain the residuals found in fitting 
a circular orbit to the radial-velocity data. The non-radial pulsation 
may be due to tidal interaction of the low-mass companion with the 
M giant. Such interaction would also cause rapid circularisation of 
the orbit for the low-mass companion. 

It is not clear what the timescale for the circularisation of the 
orbit and the dissipation of the non-radial pulsation would be. If 
the dissipation of the non-radial pulsation by convective turbulence 
is efficient and the timescale for circularisation short it is hard to 
explain why at least 25% of stars in globular clusters show LSP 
variations. We suggest that in globular clusters interactions with 
low-mass companions could be more frequent than expected. LSP 
variations then would be caused by non-radial pulsations excited 
by orbital capture of a companion or circularisation of an elliptical 
orbit. 



5 CONCLUSIONS 

We have presented simple models in order to explain the asymme- 
tries detected through infrared interferometry in the S-type symbi- 
otic star CH Cyg. We do not detect significant change of angular 
size (8.74±0.02 mas) for the M giant over a 3 year period, render- 
ing radial pulsation a less liltely explanation for the 2.1-yr variabil- 
ity in radial velocity data. We find a spherical hot dust-shell with 
an emission size of 2.2±0. 1 D, FWHM around the M giant star 
which could be responsible for some of the reported short-period 
eclipses. We find correlation between the 2. 1-yr variability and the 
variation in our closure phase. We model the closure phase as a 
large change in brightness across the M giant and/or a low mass 
companion in close orbit around the star. While the most likely 
explanation for the change in brightness can be a non-radial pulsa- 
tion we argue that a low-mass companion in close orbit could be 
the physical cause of the pulsation. The combined effect of pulsa- 
tion and low-mass companion in close orbit to the M giant could 
explain the behaviour of the radial-velocity curves and the asym- 
metries detected in the closure-phase data. If CH Cyg is a typical 
long secondary period variable then LSP variations could be ex- 
plained by the effect of an orbiting low-mass companion on the 
primary star. 
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